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High intensity acoustic radiation occurs when turbulence convects past airframe trailing 
edges. A mathematical model is developed to predict this acoustic radiation. The model 
is dependent on the local flow and turbulent statistics above the trailing edge of the flight 
vehicle airframe. These quantities are dependent on the jet and flight vehicle Mach numbers 
and jet temperature. A term in the model approximates the turbulent statistics of single- 
stream heated jet flows and is developed based upon measurement. The developed model is 
valid for a wide range of jet Mach numbers, jet temperature ratios, and flight vehicle Mach 
numbers. The model predicts traditional trailing edge noise if the jet is not interacting with 
the airframe. Predictions of mean-flow quantities and the cross-spectrum of static pressure 
near the airframe trailing edge are compared with measurement. Finally, predictions of 
acoustic intensity are compared with measurement and the model is shown to accurately 
capture the phenomenon. 
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Introduction 


Jet engines are the predominant means of propulsion for flight vehicles. These vehicles include but are 
not limited to rockets, missiles, and military and civilian aircraft. First principle mathematical models 
and physical understanding of the jet and its associated noise has eluded investigators in the fields of fluid 
dynamics and aeroacoustics. The noise from jets is problematic and linked to flight vehicle failure, hearing 
loss, community annoyance, and psychological effects. When a jet engine is integrated with the airframe of a 
vehicle, the jet flow-field can be disturbed by the airframe surface. This disturbance induces a large amount 
of additional noise production relative to an isolated jet. Accurate prediction of the noise intensity and 
associated frequencies of the jet airframe interaction will yield physical insight into the unsteady aerodynamic 
source mechanism. It will also yield a valuable tool for aircraft designers that will enable them to minimize 
this noise source. In this paper, a model is created and validated for the prediction of the power spectral 
density resulting from the jet airframe trailing edge interaction. The model is valid for a wide range of jet 
Mach numbers, temperatures, and flight vehicle speeds. 

This work draws heavily on previous experimental investigations and modeling efforts. Curie examined 
the effects of noise from surfaces by extending Lighthill’s 2 theory. Later, Ffowcs Williams and Hawkings ’ 
constructed generalized forms of the continuity and momentum equations involving two regions of the flow- 
field and created an acoustic analogy that involved the generalized density perturbation, a generalized 
Lighthill stress tensor, and two additional source terms that govern the sound from turbulence and the effect 
of moving surfaces. The latter two terms have characteristics associated with traditional acoustic monopole 
and dipole theory. Ffowcs Williams and HalT used LighthiU’s theory to make predictions of trailing edge 
noise due to eddies near to and far from the airframe edge. Soon after the work of Ffowcs Williams et 
al, ’ multiple models appeared in quick order. These included the predictive models of Jones, 1 Chase, 1, 
and notably Tam and Yu. This later work by Tam and Yu was one of the first to experimentally and 
mathematically examine the noise from a rectangular jet flow converting past a trailing edge. Spark schlieren 
images were produced of the jet trailing edge interaction. Amiet 8,9 refined a model for the airframe trailing 
edge noise due to boundary layers. Trailing edge noise prediction models based on the pioneering work of 
Amiet were developed in experimental programs at NASA Langley by Brooks and Marcolini, " Brooks et 
al . , Brooks, Brooks et al . , and Hutcheson and Brooks. Lawrence et al. modified the method 
of Amiet' 1 and included measured static pressure statistics along the trailing edge of the airframe. They 
showed that if the correct turbulent statistics were included within their model then accurate correlations 
between the airframe and observer were possible for low-speed jets. Brown created an empirical model for 
unheated subsonic jet structure interaction noise based upon master spectra and empirically fitted polynomial 
coefficients. Most recently, Goldstein et al., applied the well-known rapid distortion theory of turbulence 
to two-dimensional sheared mean-flows, and made preliminary predictions for a two-dimensional slotted jet 
flow converting past a trailing edge. However, the theory does not capture the constructive interference 
patterns that are characteristic of this phenomenon and was only demonstrated for a single jet operating 
condition and airframe position. These theories, especially those of Tam and Yu, Amiet, S/ Brooks and 
Marcolini, and Lawrence et al. form the basis of this work. 

The spectrum of jet noise is deceiving in its apparent simplicity but is extremely complicated due to 
the nature of the source. The source consists of compressible turbulence, shock- waves (for supersonic jets), 
and in the context of this paper, interactions with airframe surfaces. Figure 1 shows one such spectrum 
measured by Bridges and Brown and various predictions that correspond to the source mechanisms of jet 
noise. Sound pressure level (SPL) is on a per Strouhal number, St, basis where St is the non-dimensional 
frequency that is dependent on the fully expanded velocity and fully expanded jet diameter. The jet operates 
at a design Mach number, M ( j = 1.50, fully-expanded Mach number, Mj = 1.29, and total temperature ratio, 
TTR = 1.00. The observer (microphone) is 100 jet diameters from the nozzle exit at an inlet angle, \F = 70°. 
The spectrum due to a free jet (static) is shown as the solid line with diamonds and its associated dominant 
aerodynamic sources are labeled. Mixing noise is present at all frequencies and occurs due to the turbulent 
mixing within the jet plume. For this particular off-design supersonic jet, it is dominant in the range of 
0 + < St 0.4 and typically peaks near St « 0.3. Due to the semi-periodic shock-cell structure within the 
jet plume and the turbulence passing through it, broadband shock-associated noise (BBSAN) occurs and 
dominates the free jet spectrum in the range St 0.4. Finally, the same shock-cell structure can support an 
aeroacoustic feedback loop within the plume that induces strong discrete (screech) tones that are observed at 
St ~ 0.4 and 0.8. The dashed-dot-dot line represents a prediction of the mixing noise by Miller, 19 which is a 
modern acoustic analogy. The dashed-dot line represents a prediction of the BBSAN by Miller and Morris 
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based on the model of Morris and Miller. The combination of these predictions result in an excellent model 
for the total jet noise over a wide range of jet Mach numbers and temperatures. 

When the jet is placed near an airframe the power spectrum changes dramatically. In this example, the 
airframe (a large flat surface) extends 20 diameters downstream ( x p /D ) from the nozzle exit and extends 2 
diameters ( y p /D ) from the nozzle centerline in the cross-stream direction. These effects can be observed in 
Fig. 1 by examining the differences between the measured jet and airframe spectrum (solid line with circles) 
of Brown and the free jet spectrum (solid line with diamonds) of Bridges and Brown. For St 2) 0.2 
the total noise intensity is much lower because the mixing noise, BBS AN, and screech tone are ‘shielded.’ 
Specifically, the sound must propagate around the airframe. A prediction for jet-airframe mixing noise is 
shown as the dashed black line. This prediction is based on the preliminary work involving the tailored vector 
Green’s function of the linearized Euler equations by Miller. A prediction for the jet-airframe BBSAN is 
not shown but can easily be conducted by the method of Miller and Morris 20 combined with the tailored 
Green’s function of Miller for the airframe. At St 0.2 a large amount of additional acoustic intensity 
occurs that is due to the turbulence of the jet being deformed by the trailing edge of the airfoil. If the jet 
is farther away from the airframe, as discussed in the similarity analysis of Smith and Miller, 2 then this 
jet-surface interaction noise will not be present. The solid black line is a prediction using the mathematical 
model and physical theory developed in this paper for jet-surface trailing edge noise. 

In this paper, a model for the noise produced by the jet turbulence converting past the airframe trailing 
edge is presented. The model is based on the work of Amiet 9 ’ and inspired by Tam and Yu and Lawrence et 
al. The model depends on the integration of the source along the trailing edge of the flight vehicle airframe 
that includes the region where the jet turbulence is interacting with the trailing edge. The model seamlessly 
predicts both the noise due to jet turbulence near the airframe trailing edge and the trailing edge noise due 
to the ambient airflow (flight stream Mach number, Moo). The model scales correctly with the jet Mach 
number, jet temperature, flight stream Mach number, and is applicable to three-dimensional flows. It can 
easily be used for more complex flows from non-circular nozzles with the use of steady RANS solutions; 
however, the model presented here is independent of any external numerically generated results. 

Development of a simple empirical mean-flow model is shown. Example predictions of the empirical 
mean-flow model and its relevance to the prediction model are discussed. Auto-spectrum of the fluctuating 
aerodynamic pressure above the airframe trailing edge is compared with measurement. Predictions are 
presented and compared to measurement for five different jet Mach numbers and a heated jet. The scaling of 
the acoustic intensity is compared with measurement as the jet centerline is moved farther from the airframe. 
Finally, the validity of the model and its limitations are discussed. 

Mathematical Model 

A model is developed for the spectral density of pressure due to turbulent jet impingement and trailing 
edge noise. The model is developed based on the work of Amiet 8 ’ 9 and inspired by Lawrence et al. A 
key assumption is that the turbulence that convects past the airframe trailing edge, due to the jet flow 
and free-stream Mach number, is statistically stationary. Unlike previous theories, the noise from the jet 
impingement and trailing edge noise is due to the pressure fluctuations above the local turbulent boundary 
layer near the trailing edge. 

Above the trailing edge of the airframe the fluctuating static pressure can be written in terms of its 
wavenumber, A p = exp[i(wf — ujyi/u c — k yi yz)\. See Amiet for the analytical development of A p from a 
sinusoidal turbulent velocity impulse. Following Amiet, 9 for the pressure above the trailing edge of the airfoil, 
the method of Schwartzschild is used to calculate a Green’s function. Two solutions are obtained, one for 
quadruples and the second for dipoles, and the latter is added to the former. A transfer function (Green’s 
function) results that is used to obtain the far-field acoustic pressure by its convolution with a source. Due 
to the stochastic nature of turbulence, we must relate the far-held acoustic pressure to the turbulent pressure 
statistics (cross-spectrum) rather than choose a deterministic approach. By using the same assumption as 
Amiet, 9 that acoustic frequency is associated with a single streamwise turbulent wavenumber, we can form a 
solution that follows the same solution process as Amiet. This results in an analogous model that is similar 
to Amiet (Eqn. 4a). The resultant spectral density, S(x,lo) is, 
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S(x,u) 


Pf(x 2 -y p ) 2 u > 2 f 
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G xx l±\C\ 2 x 2 
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where c is the speed of sound, the constant Pf = (10 x 10 4//5 )~ 1 , the vector x = x(xi,X 2 ,x 3 ) is from the 
nozzle exit to the observer, x p is the streamwise distance from the nozzle exit along the centerline axis to 
a point above the trailing edge or the chord length, and u> is the radial frequency. The choice of x p will be 
apparent in the following analysis and is dependent on the trailing edge being ‘wetted’ by the jet or ambient 
flow. y p is the shortest radial distance from the nozzle centerline to the trailing edge, ds is the differential 
of the line integral that is evaluated along the trailing edge of the airframe in the spanwise direction. Note 
that the terms in the integrand are dependent on y = y{yi, y 2 , 2 / 3 )) which is the vector from the nozzle exit 
to a point on s. Figure 2 illustrates some of these variables in the coordinate system. The remaining terms 
of Eqn. 1 are now explained. 

Flight-stream propagation effects are partially taken into account by the denominator of the integrand, 


a = [(xi - yi) 2 + (3 2 (x 2 - y 2 ) 2 + {x 3 - y 3 ) 2 ] 1/2 , (2) 

where /3 is the off-design parameter, 


/3=(1 -M^) 1/2 , (3) 

and Mqo is the free-stream (aircraft) Mach number based on the local ambient speed of sound, Coo- l± is 
the local integral scale of turbulence in the cross-stream direction X 2 ■ This quantity varies along the trailing 
edge and corresponds to the cross-stream integral scale of turbulence just above the local boundary layer of 
the airfoil or within the jet plume above the trailing edge of the airfoil. We choose a model for lx as, 


i±(y,u) 



( 4 ) 


where u CyP is the local convection velocity of the coherent turbulent structures above the local boundary layer 
thickness or within the jet plume above the trailing edge. A subscript p denotes the quantity just above 
the local boundary layer thickness, <5, where S is a function of X 2 and the local flow properties. This latter 
quantity is ideally found from experiment or numerical simulation. This model is similar to that of Amiet J 
which varies based on the inverse of oj and the data of Willmarth and Roos. Here Eqn. 4 is based upon 
measured data of Brown and Wernet. 28 We propose a model for <5, 


^ _ | 0M7x p Re^ 1/o if u(y, S ) « u ^ 

[Dj (1 + aiX p /Dj) if w(y, 5) > Moo, 

where a\ = 0.07 is the spreading rate of the flow, Rei = PooUooCi/ Poo is the Reynolds number based on the 
local chord length, q, at 1 / 2 , P is the density, and p^ is the ambient kinematic viscosity. Dj is the fully 
expanded jet diameter defined by Tam, - ’ that represents the equivalent jet diameter for fully-expanded jet 
flow, 
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where D is the nozzle diameter and 7 is the ratio of specific heats. Dj is equal to D for subsonic jets. If a 
dual-stream or non-axisymmetric nozzle is present then the bypass or equivalent diameter can be used. The 
selection of Dj is important because it is related linearly to the cross-stream length scale of the turbulence 
within the jet plume above the trailing edge of the airfoil as shown in Eqn. 5. 

The auto-spectrum of the fluctuating static pressure at the surface is, 


G xx {y,u) = Sp 2 Up exp 


_ f ( 4 / 5 ) mj / i 5 + MjUj/S — ix 
V \/5 2 xl 0 4 


1 1 

- + - tanh 
8 8 


w 

20tt 


500 DjfMj, 


( 7 ) 


4 of 17 


American Institute of Aeronautics and Astronautics 


which represents an empirical fit to measured data from the experiment of Brown and Wernet. The 
functional form is very different than that used by Amiet based upon data of Willmarth and Roos. 
Alternatively, data measured directly from experiment or found through numerical simulation can be used 
in place of G xx . 

The final factor of Eqn. 1, the integration of the response within the streamwise integral of the surface 
loading (analogous to Amiet'') is, 


|£| = e 
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which has been modified from the form originally proposed by Amiet. /i is, 
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The function E* of Eqn. 8 is an integral of the form, 
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Equation 11 is in a well-known form that can easily be rewritten in terms of Fresnel integrals, 
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The Fresnel integrals can be written as a summation of two Error functions 
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The argument of the error function appearing in the expansion of F s and F c is complex. The error 
function is available in tabular form and can easily be evaluated with many modern computing languages. It 
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is often convenient (or necessary) to evaluate an error function with a complex argument using the Faddeeva 
function, 


w(z) = exp [— z 2 ] erfc[— iz], (18) 

where erfc(,z) is the complementary error function with complex argument 2 . Evaluation of the Faddeeva 
function can be performed numerically with the algorithm developed by Poppe and Wijers. erf is related 
to w(z) by, 


erf(z) = w[—iz] exp[^ 2 ] — 1. (19) 

Equation 1 and its arguments are defined except for the mean-flow quantities. These quantities can easily 
be found with experimental or numerical investigations. Here, we will propose a simple empirical model for 
single-stream jets. The mean-flow model is valid for jet structure interactions in the range of 0 + < M :} <1.4 
and 1 < TTR ;$ 3.2 and should be used with caution otherwise. It is in-part based on the models of Lau et 
al. and Lau. The airframe structure has some effect on the jet mean-flow and development but for now 
we will assume that the plume development is axisymmetric. The core length x c of the jet is, 


x c = Dj {4.2 + — Moo) 2 } , 

and the streamwise velocity component on the jet centerline axis is, 
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The streamwise velocity component relative to the nozzle exit is, 
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where r = \J y\ + y\ and rf = (r — Dj /2) /y\ . The error functions have real arguments and the use of the 
Faddeeva function relation of Eqn. 18 is not required. The spreading rate of the jet is, 


Sr, = 0.177 { [l - 0.294M 2 ] [l + 0.50 (M 2 - 1)(T,-/T 0 - 1.40) 2 ] } , (24) 

where T a is the stagnation temperature in the plenum of the nozzle. The temperature is required to find the 
Mach number. A similar equation for T is created for this purpose, 
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where S v t = 1A55 V and a t = y/Tr/8 rit is the Gortler spreading parameter. Quantities important for the 
developed model are the convection speed, 
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the temperature T p = T(x p , \Jy p + 2 / 3 ), local speed of sound c p = j/^/RXp, and Mach number, M p = u p /c p . 
The convection speed of the turbulence is assumed to be approximately u CjP = 0.70it p with associated axial 
wavenumber k x = oj/u CtP . The later assumption is required in the formation of Eqn. 1. 
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Results 


This section discusses the coordinate system in relation to the previously defined independent variables 
within the mathematical development section. An example mean-flow prediction is shown based on the 
empirical model and compared with a steady RANS solution generated by a separate numerical method. 
The auto-spectrum model is compared with measurement. Then, example predictions are shown for a number 
of jet operating conditions and plate locations. 

Figure 2 illustrates the problem geometry and coordinate system. The origin of the coordinate system 
is the nozzle exit. The aq-axis coincides with the centerline axis of the nozzle, the ^ 2 -axis is normal to and 
points away from the airframe, and the a.' 3 -axis is in the cross-stream direction. M 0 0 is in the direction X\. 
y p is the negative distance from the centerline axis to the trailing edge bounded on the aq-axis and x p is the 
distance from the nozzle exit to the trailing edge bounded on the aq-axis. The source vector, y, is defined 
from the nozzle exit to a position along the airframe trailing edge defined by the line s. The differential 
element, ds, of the model equation is bounded on s and y is smallest at y = y(x p , —y p , 0). Note the trailing 
edge of the airframe is not restricted to a constant y i position and can vary in space, x is a vector from the 
nozzle exit to the observer. T is the angle between the downstream centerline nozzle axis and the vector a:, 
and <j> represents the azimuthal angle about aq. Observers with angle 0 = 0 are restricted to the aq — aq 
plane. 

Tables 1 and 2 summarize the positions of the airframe and jet operating conditions, respectively, that 
are available to the author for validation. For each jet operating condition (row of Table 2), Table 1 shows 
various airframe positions relative to the nozzle exit. For every streamwise plate position, x p /D, multiple 
cross-stream plate positions, y p /D, are available. The first column of Table 2 states the nozzle name, Small 
Metal Chevron (SMC) and a three digit series number. Note neither nozzle has chevrons. Nozzle pressure 
ratio (NPR) is the ratio between the total pressure in the nozzle plenum and the ambient static pressure. 
The other quantities of Table 2 are defined in the previous section. 

An example calculation of the variation of the mean-flow along the trailing edge, s, is shown in Fig. 3. 
The mean-flow quantities are induced by a Mj = 0.50 and TTR = 1.00 jet flow, x p /D = 10 and y p /D = —2 
airframe position, and M x = 0. This condition corresponds to the first row of Table 2. The lines represent 
predictions using the developed empirical mean-flow model for single-stream jets. The symbols are from 
a respective steady RANS solution that includes the flat plate geometry and SMC000 nozzle relative to 
the plate, u, T, and M are the streamwise velocity component, static temperature, and Mach number 
respectively above the local boundary layer thickness or above the local jet length scale in the X2 direction. 
The solid line and diamonds represent u and are predicted from the empirical model (Eqn. 23) and the 
steady RANS solution respectively. For this particular case the induced velocity due to the jet, according to 
the empirical model and the steady RANS solution, is less than 1 m/s. The largest magnitude of u agrees 
very well with the steady RANS solution and both form a nearly Gaussian shape. The dash-dot line and 
circles represent T and are predicted from the empirical model (Eqn. 25) and the steady RANS solution 
respectively. The ambient temperature is 293 K. Because the jet is unheated, T is less than the ambient 
temperature within the jet flow. The largest change in T relative to the ambient is 0.1 K according to the 
model and 0.5 K according to the steady RANS solution. Like u, the shape is Gaussian but inverted about 
the ambient temperature. Heated or high speed cases induce very large variation in T such as Mj = 1.50 or 
Mj = 0.68 (TTR = 1.93) in Table 2. The final mean-flow comparison shown in Fig. 3 is for M. The dash- 
dot-dot line and triangles represent M and are predicted from the empirical model and the steady RANS 
predictions respectively. M is calculated by the ratio of Eqn. 23 and the speed of sound based upon Eqn. 25. 
In this case, the agreement between the model and steady RANS solution has the same characteristics as 
the agreement in the u comparison. 

Outside the range of \x^/D\ ;$ 2, u, T, and M are at nearly the ambient flight stream condition. Within 
the disturbed region of the mean-flow, the dominant noise is due to the jet structure interaction and outside 
this region traditional trailing edge noise dominates. These quantities and others continuously transition and 
the resultant sources of the jet-surface interaction and traditional trailing edge noise change with varying 
jet operating conditions, airframe position, and free-stream Mach number. Note that the quantities u, T, 
and M can easily be provided by a steady RANS solver or another approach for potentially more accurate 
predictions. The induced quantities on the airframe trailing edge can be very small as illustrated in this 
case. Thus, the simple symmetric model is useful for many calculations and it will be shown that it works 
well for cases where the jet is highly deformed. One very important point must be stressed: the convection 
velocity and convective Mach number of the turbulent structures of the jet are much larger than the induced 
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velocity and Mach number at the trailing edge. 

The variation of the auto-spectrum of the fluctuating pressure above <5, G xx , as predicted by Eqn. 7 is 
compared with measurements in Fig. 4. The jet operates at Mj = 0.50 and TTR = 1.00 and the ambient 
Mach number is = 0.001. The airframe position is y p /D = —2 and the streamwise trailing edge location 
Xp/D is located at 4, 8, 10, and 15. 

The measurements were conducted by Brown and Wernet and the data were processed by Palumbo. 
Measurements at other conditions shown in Table 2 are not available to the author for validation of G xx . 
The pressure sensor was located as close to the trailing edge as possible without compromising the structural 
integrity of the plate. For more information on the measurements see Brown and Wernet. 

Predictions using Eqn. 7 are shown as a solid, solid-dot, solid-dot-dot, and dashed lines for x p /D = 
4, 8, 10, and 15 respectively. Processed measurements are shown as diamond, circle, triangle, and square 
symbols at the same respective locations as the predictions. There is a large amount of scatter present in the 
measured data. Measurements show that as x p /D increases, the magnitude and the spectral width increase. 
The frequency of the peak magnitude decreases with increasing x p /D. Lower frequency spectral content 
increases in magnitude with increasing x p /D and is due to the growth of the turbulent scales with increasing 
streamwise distance. At large downstream locations, such as x p /D = 15, the overall width increases but 
intensity decreases. 

These trends are generally captured by the proposed model, Eqn. 7. The model is linked to the jet 
operating condition and airframe position through a semi-empirical model of the statistical flow quantities. 
Predictions contain a low-frequency fall-off. The predicted high-frequency fall-off is too large relative to 
measurement, and a log-normal model can correctly capture this fall-off. The form of G xx is intended not 
to capture the log-normal fall-off of the measurement because the model must reflect the energy of noise 
radiated from a distance S above the airframe trailing edge, and not the energy on the wall. The scaling 
terms present in G xx for Mach number and temperature effects are based purely on theory. It will be shown 
in the following noise predictions and comparisons that the choice of scaling terms is valid for the proposed 
model of G xx . 

Jet structure interaction noise predictions using Eqn. 1 are now compared with measurements of Brown. 
Predictions have been performed for all the airframe positions and jet conditions shown in Tables 1 and 2 
respectively. The first prediction is shown in Fig. 5. Predictions of SPL per unit St are compared with 
measurements from the Mj = 0.50 and TTR = 1.00 jet. The nozzle is convergent M,i = 1.00 (SMC000) and 
the plate location is x p /D = 10 and y p /D = —1. Recall that the angle, T, is measured from the upstream 
axis and is shown in each sub-figure. Free jet measurements of Bridges and Brown are shown as a line with 
diamond symbols and are labeled ‘SHJAR’ (Small Hot Jet Acoustic Rig). Jet with airframe measurements of 
Brown are shown as a line with circles and are labeled l JSIT’ (Jet Structure Interaction Test). Predictions 
of the model equation are shown as a solid line without symbols. The free jet measurements (SHJAR) are 
shown along with the jet with airframe measurements (JSIT) to illustrate the spectral changes that occur 
when placing an airframe surface near the jet plume. Specifically, the free jet mixing and shock-associated 
noise can change dramatically by the presence of the airframe surface. Also, additional sources associated 
with the jet structure interaction noise are created by the introduction of the airframe into the flow. The 
predictions illustrate one type of jet structure interaction noise that is due to the turbulence being deformed 
by the plate trailing edge and its subsequent propagation via incident and diffracted waves. Predictions in 
this paper, particularly at low through mid St, should be compared with the jet with airframe measurements 
(JSIT). The presented results are corrected for atmospheric absorption (lossless spectra) and corrected to 
a distance of 100 nozzle exit diameters ( R/D = 100). Note that correcting for the atmospheric losses and 
propagation distance can be difficult for jets near airframes because the source is highly distributed and the 
distance that the acoustic waves travel to the observer can be difficult to calculate, especially so when the 
propagation mechanism primarily includes diffraction. 

At the various observer angles the predicted spectral magnitude, spectral width, and fall-off at high 
and low frequencies is satisfactory predicted. However, in the upstream direction the spectral magnitude is 
lower than measurement by 4 dB. The low-frequency fall-off matches measurement excellently in all cases, 
and this is observed even below St ~ 0.01 when data are available. High frequency fall-off is also very 
satisfactory and there is little to no noise contribution above St ~ 0.8. This is highly desirable because 
the spectrum above St « 0.8 is dominated by the jet mixing noise, as shown and discussed in the previous 
section. In the sideline and upstream direction multiple peaks can be observed both in the prediction and 
measurement. These peaks are due to the diffraction of the noise by the trailing edge of the airframe and 
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are taken into account by Eqn. 8. The first and dominant peak frequency, for this case, generally matches 
that of measurement. As the observer angle is increased the peak magnitude increases in frequency and the 
multiple broad spectral peaks are unified. 

A second comparison is shown in Fig. 6 where Mj = 0.68 and TTR = 1.93. The nozzle is convergent 
Md = 1.00 (SMC000) and the plate location is x p /D = 20 and y p /D = — 1. Relative to the previously 
presented case the jet Mach number and temperature has increased and the plate position is extended. 
Multiple broad spectral humps are again observed both in the predictions and measurement and are more 
prevalent in the upstream and sideline observer directions. Jet mixing noise is more dominant at lower 
frequencies (St ~ 0.5) as the plate is extended. Predictions at angles 4' = 110 through 150 deg. match 
measurement. At 4 f = 50 deg. the prediction is approximately 8 dB lower than measurement in peak 
magnitude and lower in frequency by St ~ 0.03. However, in this direction the multiple spectral peaks 
due to diffraction by the trailing edge are most apparent both in prediction and measurement. Although 
the peaks are not coincident in frequency with measurement, the respective predicted and measured peaks 
are present. Trends in prediction and measurement show that the narrower spectral width and increased 
magnitude in the upstream direction is due to a higher value of x p /D and not due to jet Mach number or 
increased temperature. 

Figure 7 shows predictions and measurements of the Mj = 0.98 and TTR = 1.00 jet. The nozzle is 
convergent and the plate location is at x p /D = 4 and y p / D = —2. Relative to the previous cases the nozzle 
is located very close to the plate edge. Also, comparing the free jet (SHJAR) and jet and airframe (JSIT) 
measured spectra relative to previous cases, shows that in the sideline through upstream direction there is 
very little change. At low- to mid-frequencies in the upstream and sideline directions there is only a very 
small increase in measured acoustic intensity due to the jet structure interaction. The predictions at these 
angles agree extremely well with measurement and correctly predict this small increase. The difference at 
higher frequencies is due to a type of shielding that is correctly predicted by Miller and Smith. At the 
downstream angles, 4' = 130 and 150 deg., the jet mixing noise dominates any contribution from the jet 
structure interaction noise. As expected, the predictions are completely dominated by the jet mixing noise. 
By comparing the predicted spectral shapes relative to previous cases shows that there are still multiple 
broad lobes present, however, the second broad loab is dominant in magnitude. 

The next comparison is shown in Fig. 8 for an off-design Mj = 1.29 and TTR = 1.00 jet. This case is 
very similar to the one in the introduction as shown in Fig. 1. The nozzle is convergent-divergent Md = 1.50 
(SMC016) and the plate location is at x p /D = 10 and y p /D = —2. Free jet spectra show all the characteristics 
of coherent and incoherent mixing noise and BBSAN. With the addition of the plate the measured spectra 
changes dramatically. These sources and changes can be predicted as discussed in the introductory section. 
Overall, mixing noise and BBSAN dominate all jet structure interaction noise above St « 0.4. At the 
upstream and sideline direction, the peak jet structure interaction noise is predicted correctly in terms of 
magnitude and peak frequency. At lower frequencies and at T = 50 deg., the broad lobe is predicted 
correctly at St ss 0.02 but is 4 dB too low. The total noise at low frequencies in the downstream direction 
is a combination of the coherent jet mixing noise and jet structure interaction noise. 

A final comparison of predictions relative to varying observer location is shown in Fig. 9 for the Mj = 1.50 
and TTR = 1.00 jet. The nozzle is convergent-divergent Md = 1.50 (SMC016) and the plate location is at 
Xp/D = 20 and y p /D = — 2 . The jet operates supersonically and on-design. Unlike the previous case, the 
measured spectra only shows extremely weak BBSAN near St ss 0.8. At the sideline location, the predicted 
interference pattern is predicted but is 0 to 6 dB under-predicted (at the ‘troughs’). Unfortunately, at 
upstream angles the under-prediction is worse. At 4' > 110 deg. the magnitudes and peak frequencies are 
predicted more accurately relative to the upstream angles. For example, at the extreme downstream angle 
the peak frequency and fall-off at high and low- frequencies is predicted correctly, relative to measurement. 
These previous two cases presented illustrate the on-set of compressible turbulence. This on-set has been 
mainly incorporated in this prediction method through empirical correction. Excellent measurements of 
the statistics of compressible turbulence and its onset for high speed jets near airframe surfaces can guide 
theorists to create superior predictions. 

We will now examine the variation of the jet structure interaction noise by varying the radial position of 
the airframe. Figure 10 shows Mj = 0.68 and TTR = 1.93 jet spectra at 4< = 90 deg. and R/D = 100 from 
the convergent SMC000 nozzle. The plate streamwise location is x p /D = 10 and y p /D is varied from -1 to 
-6. Each set of spectra consist of the free jet and JSIT measurement and associated prediction. Note that 
the free jet spectrum is the same for each plate position. The plate locations are labeled on the left and the 
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maximum predicted SPL per unit St is shown above for each set. 

As the plate is moved farther from the jet centerline axis the jet structure interaction noise intensity 
decreases. This can be observed by comparing the JSIT and SHJAR measurements in Fig. 10. Predictions 
at each y p /D location show the same relative trend as the SHJAR measurement. Also, the predictions at 
each y p /D location compare favorably relative to the JSIT measurement and illustrate decreasing amplitude 
and St with increasing jet structure separation. 

The final comparisons of this paper are shown in Fig. 11 for the Mj = 1.29 and TTR = 1.00 jet. These 
comparisons are like the previous except the convergent-divergent M,i = 1.50 (SMC016) nozzle is used 
and the jet is off-design and unheated. Examining the relative difference between the JSIT and SHJAR 
measurements shows similar trends to the heated subsonic case except for the strong presence of BBSAN. 
The low- and high-frequency fall-off and peak magnitudes show the same relatively good agreement with 
the JSIT measurements as previous cases. In this particular case, there are strong constructive peaks that 
are not observed in the measured spectrum, and this can possibly be attributed to data processing. The 
predictions relative to measurement in Figs. 10 and 11 strongly suggest that the model is capable of correctly 
predicting the fall-off of intensity and peak frequencies for subsonic, supersonic, and cold and hot jets. In 
summary, the model is validated for single-stream jet flows that operate from Mj = 0.50 through Mj = 1.50 
and TTR = 1.00 through TTR = 1.93. There is no reason to expect that the intensity scaling will be 
incorrect for lower jet Mach numbers or higher temperatures. 

Conclusion 

A mathematical model is proposed for the prediction of noise from jet structure interaction. The model 
can seamlessly predict the jet structure interaction noise and noise induced by ambient air passing the trailing 
edge. A simple mean-flow model is developed and validated with steady RANS solutions, but steady RANS 
solutions based on more realistic nozzles can also be used for more accurate predictions. Furthermore, a 
model for the auto-spectrum of fluctuating pressure above the boundary layer is proposed and validated with 
measurement. The auto-spectrum model is dependent on the jet operating condition and local mean-flow. 
The jet structure interaction noise model is validated successfully for a wide range of available single-stream 
jet Mach numbers, temperature ratios, and airframe positions. Predictions compare very favorably with 
measurements at various observer angles. The decay of acoustic intensity with increasing jet centerline 
airframe distance also compares favorably. 

The model is expected to scale correctly outside its range of validation with Mj and TTR and can 
easily be extended to more complicated jet flows. However, careful measurements of compressible turbulent 
statistics within heated high-speed jet plumes are necessary to make successful predictions at higher jet Mach 
numbers and temperatures. These measurements will aid in the creation of a more general cross-spectrum 
model of fluctuating pressure. The author is hopeful that in the future such difficult measurements or even 
experiments will be undertaken. Finally, the combination of this source model with other models of the 
author yields an almost complete prediction method for off-design, heated, compressible turbulent jets. 
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Tables 


Table 1. Plate Locations. 


Axial Distance ( x p /D ) 

Radial Distance ( y p /D ) 

4 

-1, -2, -4, -6, -8, and -10 

10 

-1, -2, -4, -6, -8, and -10 

20 

-1, -2, -4, -6, -8, and -10 


Table 2. Jet Operating Conditions. 


Nozzle 

M d 

NPR 

TTR 

Uj (m/s) 

Tj (K) 

Mj 

SMC000 

1.00 

1.186 

1.00 

167.4 

279.0 

0.50 

SMC000 

1.00 

1.363 

1.93 

310.1 

517.6 

0.68 

SMC000 

1.00 

1.850 

1.00 

312.9 

244.3 

0.98 

SMC016 

1.50 

2.733 

1.00 

383.4 

219.8 

1.29 

SMC016 

1.50 

3.671 

1.00 

427.4 

202.1 

1.50 
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Figure 1. Predictions and measurements of the Md = 1.50, Mj = 1.29, TTR = 1.00 jet at a distance of 100 jet diameters, 
i/> = 70 deg., on the opposite side of the airframe from the jet. The jet-airframe location is at x p /D — 20 and y p /D = —2. 
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Figure 2. The coordinate system illustrating the nozzle off-set from the trailing edge of the airframe. 
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Figure 3 . Predicted (lines) mean-flow quantities induced along the trailing edge from a Mj — 0.50 and TTR = 1.00 
jet flow, Xp/D — 10 and y p /D — —2 airframe position, and M m = 0. The predictions are compared with steady RANS 
solutions (symbols) that include the flat plate geometry. 
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Figure 4. Predictions of G xx with measurement at the trailing edge induced by a Mj = 0.50 and TTR = 1.00 jet flow 
with varying x p /D and y p /D = —2 airframe position, and = 0.001. 



Figure 5. Predictions of SPL per unit St compared with measurements from the Mj = 0.50 and TTR = 1.00 jet. The 
nozzle is convergent Mj = 1.00 (SMC000), the plate location is x p /D = 10 and y p /D = — 1, and microphones are at 
R/D = 100. 
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Figure 6. Predictions of SPL per unit St compared with measurements from the Mj = 0.68 and TTR = 1.93 jet. The 
nozzle is convergent Md = 1.00 (SMC000), the plate location is x p /D = 20 and y p /D = — 1, and microphones are at 
R/D = 100. 



Figure 7. Predictions of SPL per unit St compared with measurements from the Mj = 0.98 and TTR, = 1.00 jet. The 
nozzle is convergent Md = 1.00 (SMC000), the plate location is x p /D = 4 and y p /D — —2, and microphones are at 
R/D = 100. 
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Figure 8. Predictions of SPL per unit St compared with measurements from the Mj = 1.29 and TTR = 1.00 jet. The 
nozzle is convergent-divergent Md = 1.50 (SMC016), the plate location is x p /D = 10 and y p /D = — 2 , and microphones 
are at R/D = 100 . 



Figure 9. Predictions of SPL per unit St compared with measurements from the Mj = 1.50 and TTR, = 1.00 jet. The 
nozzle is convergent-divergent Md — 1.50 (SMC016), the plate location is x p /D = 20 and y p /D = — 2 , and microphones 
are at R/D = 100 . 
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Figure 10. Predictions of SPL per unit St compared with measurements from the Mj = 0.68 and TTR = 1.93 jet. The 
nozzle is convergent Md = 1.00 (SMC000), the plate location is x p /D — 10 and y p /D is varied, and the microphone is at 
R/ D — 100 and = 90 deg. 



Figure 11. Predictions of SPL per unit St compared with measurements from the Mj = 1.29 and TTR = 1.00 jet. 
The nozzle is convergent-divergent Md = 1.50 (SMC016), the plate location is x p /D = 10 and y p /D is varied, and the 
microphone is at R/D — 100 and ^ = 90 deg. 
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